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Transverse parton distribution functions at next-to-next-to-leading order: the 

quark-to-quark case 
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We present a calculation of the perturbative quark-to-quark transverse parton distribution func- 
tion at next-to-next-to-leading order based on a gauge invariant operator definition. We demonstrate 
for the first time that such a definition works beyond the first non-trivial order. We extract from 
our calculation the coefficient functions relevant for a next-to-next-to-next-to-leading logarithmic 
Qt resummation in a large class of processes at hadron colliders. 
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PACS numbers: 12.38Bx 

Parton distribution functions (PDFs) describe the mo- 
mentum distribution of quarks and gluons inside hadrons, 
and are essential inputs for physics program at hadron 
colliders such as the CERN Large Hadron Collider 
(LHC). Their usefulness resides in the factorization theo- 
rems [1| which separate physics at different length scales. 
While PDFs may depend on the 4-momentum of the par- 
ton (the so-called fully-unintegrated PDFs), most physi- 
cal observables are only sensitive to the forward compo- 
nent of the parton momentum. To be more precise, it 
is convenient to introduce two light-like vectors n and n 



satisfying n 



and n-n = 2, where n is along the 



beam direction of the hadron. Any 4- vector v 11 can then 
be decomposed as = h ■ v /2 + n ■ v + v±. For 
an energetic parton along the n-direction, its momentum 
q 11 has the hierarchical structure q = (n ■ q,n ■ q,q±_) ~ 
Q(1,A 2 ,A), with A <C 1. For hard interactions, the 
smaller components n ■ q and q± can often be neglected. 
The resulting PDFs are referred to as "collinear PDFs" , 
which only depend on the fraction z — n ■ q/n ■ p with p 
being the hadron momentum. 

The collinear PDFs can be defined rigorously in quan- 
tum chromodynamics (QCD) as matrix elements of cer- 
tain non-local operators 2]. For example, the quark PDF 
is given by 

x (N(p)\$W](tn) | [WH](0)\N(p)) , (1) 

where ip is the quark field and W is a light-like Wilson line 
which renders the operator gauge-invariant. The PDFs 
are non-perturbative objects, but their renormalization 
group (RG) evolution can be determined perturbatively 
as long as the factorization scale /i ;§> Aqcd- Their 
anomalous dimension functions (splitting functions) have 
been computed up to 3 loops 

Certain observables, such as the transverse momen- 
tum distributions in the production of the Higgs bo- 
son or gauge bosons, however, are sensitive to the 
perpendicular component q±_ of the parton momen- 
tum. One therefore needs to keep q± in the defini- 
tion of the relevant PDFs, which are referred to as 



transverse-momentum-dependent PDFs (TMDPDFs) or 
simply transverse PDFs. Closely related to the trans- 
verse PDFs is the Collins-Soper-Sterman (CSS) formal- 
ism of transverse momentum resummation [ij , which ad- 
dresses the divergent behavior of fixed-order calculations 
at small transverse momentum. The CSS formula for the 
Drell-Yan process can be written in the form 



da 



dQ 2 dQ 2 T dy 9Q 2 s J (2tt 



d 2 x± 




A(a s {p,)) +B(a s (fl)) 



V q / Nl (zi,x T ) Vq/ N2 (z 2 , x T ) + (q O q) 



(2) 



where is the momentum of the Drell-Yan pair, x± 
is the variable conjugate to Q± (also referred to as the 
"impact parameter" b in the literature), and /if, = bo/xr 
with bo = 2e~ 7E . The functions V can be interpreted as 



transverse PDFs. 



For x\ 



-x 2 , < 



match V onto the collinear PDFs 



V q / N (z,x T ) 



*(/*&)) fc/Niz/ZiPb) 



with perturbatively calculable coefficient functions C q {. 
These coefficient functions were determined to the next- 
to- leading order (NLO) in The general structure of 

at NLO was obtained in [6] including C g t relevant for 
Higgs production. Previous results beyond NLO were 
based on a modification of the CSS formula 171, intro- 



ducing functions H. 



ijl—ab ? 



which are related to the con- 



volutions of Ci a and Cjb functions. The next-to-next-to- 
leading order (NNLO) corrections to ri were computed 
[8| for the Drell-Yan process and for Higgs production. 
The method used in these results is to assume the formula 
©, and compare its fixed-order expansion with explicit 
calculations of the Qt spectrum. It is however desirable 
to have an operator definition of the transverse PDFs, 
and compute the matching coefficient functions from the 
definition. 

One may extend the definition ([I]) to the case of trans- 



2 



verse PDFs as [9] 



dt 
2^ 



-iztn-p 



x (N(p)\[i>W](tn + x ± )^[W^](0)\N(p)) . (3) 

This definition, however, is problematic since new diver- 
gences associated with the light-cone propagators arise, 
which are not regulated in dimensional regularization. 
Therefore, one needs to supplement Eq. ([3]) with some 
extra regulator to make it well-defined. Generically, the 
result will then depend on this extra regulator, while in 
physical observables such dependence will cancel. 

Different kinds of regulators have been proposed in 
the literature. The original paper Q employed a non- 
light-like axial gauge. That was employed in subsequent 
NLO calculations [f0| • More recently, Collins introduced 
a gauge-invariant definition utilizing non-light-like Wil- 
son lines If L An equivalent [l2| definition has been put 



forward in [13J] . Variants of Smirnov's analytic regulator 



were used in [t| [f5|, [f6| • All these approaches were 
argued to be valid to all orders in perturbation theory 
based on factorization properties. However, explicit cal- 
culations were only carried out at NLO. In this Letter, 
we report the first NNLO result for the transverse PDF 
in the quark-to-quark case, using the regulator proposed 
We also extract the NNLO coefficient function 



15 



Cqq , which is the first direct calculation of this function. 

We consider processes where a qq pair annihilates 
into some color neutral final state F: q(pi) + q(p 2 ) 
F(Q), with pi along the n-direction and p 2 along the in- 
direction. We follow closely the formalism in [9J, where 
the transverse PDF for quarks along the n-direction is 
defined as in Eq. © , while the PDF for quarks along the 
opposite direction Bq/M is defined with n f-> n. To com- 
pute the matching functions, we replace the hadron field 
N with a quark field and evaluate the matrix element 
in Eq. ([3]). The B functions require an extra regulator 
beyond dimensional re gula rization, for which we adopt 



the one introduced in |15| . namely, we multiply a fac- 



tor (v/n ■ k) a for each emitted parton with momentum 
k. The B functions will contain poles in the analytic 
regulator a, which however will cancel in the product 
B q /q Bq/q. A remnant of this regulator dependence is the 
collinear anomaly resulting in a dependence on the 
hard momentum transfer Q 2 in the product, which can 
be refactorized as: 



[Bq/q(zi, X%, /i) Bq/q(z 2 , Xt, h)] 



(4) 



„2 n 2 \ -Fq<i(L±4i) 

Jbrp \c£ 



4 e -27E 



Bq/q(z 1 ,L±, Bq/q(z 2 ,Lj_, (J,) 



where L± = ln(x|,/i 2 /&Q). Note that B q / q = Bq/q, we 
therefore do not distinguish them anymore. The Bi/j 
functions can be regarded as the process-independent 



transverse PDFs, which can be matched onto the 
collinear ones for xt <C 1/ Aqcd in the form 



Bi/j{z,L A 



.a0 = £ 



(•5) 



We define the perturbative expansion of the Ii/j func- 
tions as Ii/j = J2n( a s / (^ :7r )) n ^ijj > an d similarly for other 

functions. The second order coefficient / / is the main 

q ' q . 

result in this Letter, while the other combinations of i 
and j will be presented in a forthcoming article. 

At leading order (LO), it is clear that B q °^(z, x\, /i) = 

B^-{z, x\, /i) = 5(1 — z). At NLO, these functions were 
calculated in [t| using a slightly different way of regu- 
larization. For our purpose, we need to recompute them 
with the regularization scheme of f5[ . The results read 



B^(z,x 2 T ,fi) = c F e^ L ^ +2a ^ E 



ill 

n-pi 
A 4 

1(1) /_ „2 



r(-e - a) 
T(f + a) 



- (f-z)- 1+ "[4z + 2(f- £ )(f-z) 2 ], 
B£/ 9 (z,a%,VL)=C F e^-^T(-e) (6) 



n-p 2 



{l-z)- 1 - a [Az + 2(l-e)(l-zY 



It is easy to check that the poles in a vanish in the com- 
bination 5(1 - zi) B^]q(z2) + B^ q (zi) 5(1 - z 2 ). The bare 

F„^ and I, functions can then be extracted from this 
ii ill 

combination following Eqs. (j4]) and ([5]). The remaining 
poles in the dimensional regulator e = (4 — d)/2 can be 
renormalized in the MS scheme: 



F^(L ± )=Z^)+F il (L ± ,^) 



(7) 



Note that for the renormalization at NNLO, we will need 
Iq) g and Zgfq \ which we also computed. 

At NNLO, the transverse PDFs receive 3 classes of 
contributions. In Figure [I] we show a typical Feynman 
diagram for each of them: (a) virtual+real diagrams; (b) 
double gluon emission diagrams; (c) quark-antiquark pair 
emission diagrams. The virtual+real diagrams are rel- 
atively easy to calculate. After carrying out the loop 
integrals, we encounter familiar integrals which already 
appeared at NLO. The virtual+real diagrams contain di- 
vergences requiring coupling constant renormalization, 
for which we include diagrams with the loop replaced 
by a counter-term. The main complication comes from 
the double real emission diagrams (b) and (c), where two 
propagators are raised to non-integer powers. We man- 
aged to reduce them to two-fold integrals involving hy- 
pergeometric functions. We then perform the remaining 







(a) 



(b) 



(c) 



FIG. 1. Sample Feynman diagrams for our calculation: (a) virtual+real; (b) double gluon emission; (c) quark-antiquark pair 
emission. 



integration by a systematic expansion in a and e. To 
this end we made extensive use of properties of hyper- 
geometric functions and the HypExp package [l7|. The 
anti-collinear PDF B q / q can be obtained similarly, with 
n ■<-> n. Note however that the analytic regulator is still 
given by (yjn ■ k) a . 

The B q ) and B?) q functions contain poles in the an- 
alytic regulator a. It is important to verify that these 
poles cancel in the combination 6(1 — z\) B- L(z2) + 

B^ q ( Zl ) 6(1 - z 2 ) + B ( q f q ( Zl ) B^-(z 2 ). We therefore show 
below the singular structures in these two functions: 



BfUz,L A _) = AC F e 2 ^ L ^ 



2a / \ 2a 
n ■ pi \ IV 
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C, 



1 

3 + 



A 4 / VM, 
2C F (l-z) 5(1 -z) 



e a 



■■ > ! I ^ + C2 



4 
7s 



fCa 
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C A 



67 101 23CsV„ (5 28 
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2? 

(8) 



Po ( 1 



+ (2 



Bf/ q (z ) L x )=4C F e 2 ^ 



+ B q 2 / ) q ct (z,L ± ) + O(a ) 

\ 2a / N 2a 
/I \ I V 



P£\ Z ) ( 1 



^ + C2 ) + 



n-p2/ 

2C F (l-z) 6(1 -z) 
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) 



^ 2 / ) 9 ct (z,L ± ) + O(a ), 

>(0), 



(9) 



where P^(z) = 2C F [(1 + z 2 )/(l - z)}+, (3 = HC A /3 - 

3 (2),c 

<?/<? 



AT F nf/Z. The function £>( 2 ; )j ct comes from a, renormal- 



ization and is given by 



B™*\z,L L )=iC„e 



a V 2 



(10) 



From the above formulae, it is clear that the pole 
terms with the color structures C f Ca and CpTpnf van- 
ish in the sum 6(1 - Zl )B^ q (z 2 ) + B ( q J q (z l )6(l - z 2 ). 
The remaining singularities are canceled by the product 



a), 



(2) 

We are now ready to extract the functions and 

('2) 

, following the procedure of Eqs. ((4]) and ([5]) and car- 
rying out the renormalization as in Eq. ([7]). We have 

t 2 \ 

checked that the F qq ' function extracted from our cal- 
culation agrees with the expression given in [9] and that 
the / 



q l q function satisfies the RG equation 
-Li A*) 



dI q / q (z,L 



din [i 



r c Usp ( a s)i± - 2 1 q (a s ) I q/q (z,L±,Li) 



k 



2^ ^-I q/k (^L ± ,fx)P kq (z/^fi), (11) 



where the anomalous dimensions r^ IO „ and j q up to 2 



cusp 



loops can be found in [9J , and the splitting functions Py 
up to 2 loop order can be found in [18J. The finiteness 



(2) 

of I q i q and its RG properties demonstrate, for the first 
time, that the operator definition for the transverse PDFs 
supplemented with the analytic regulator is valid beyond 
the first non-trivial order. 

(2) 

We finally give the scale-independent part of the I q ^ q 
function, which is the main result of this Letter. It can 
be written as 



266C 3 



6(1 -z) 



C F — 7- + C f Ca 



/3032 67C 2 



5C4 + C F T F n f 
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TpTlf 
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68 4$, 32 4 4 
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(2 - 24z)tf + (3 + 7z)H a . Q + 2(1 + z)ff ,o,o 
2zH 1 + {l-z) (6C2 - 22 + 4ff ,i + 12ffi,o 



CfCa 



(2 + 10z)iJ - 4ziJ ,o - 2ziJi 



44 



+ (1 -*)( — - 66 - 4F lj0 



C F T F 



-50 + 38z 20 + 8z 2 - 22z 

n + n H ° + 5 #0,0 



+ 4(1 + 2)^0,0,0 



-C F T F n f (l - z) 



(12) 



where Pf g \z) = 2T F [z 2 + (1 - z)% P$(z) = 2C F [1 + 
(1 — z) 2 ]/z, and i?{ m } = H({m},z) are harmonic poly- 
logarithms introduced in [191 ] . 

The coefficient function Cm, relevant for the Drell-Yan 
process is related to I g i q by [9( 

C qq (z,a s (n b )) = \C v (-f4,fib)\lq/q(z,0,a a (ni)) , (13) 

where the function CV is determined from the virtual 
corrections to the Drell-Yan process and can be extracted 
from 20] . The NNLO expression for C qq reads 

C$(z, aM) = l ( q f q (z,0) + C 2 F (l- z) (14& - 16) 
'255 _ „„„ 87C 4 ' 



+ C F d(l-z) x 



19C 2 - 30C 3 



C A - 
T F n f 



51157 IO6IC2 313C 3 



648 18 9 
4085 182C 2 4C 3 



-80 



162 



9 



(14) 



Starting from this expression, we have checked that we 
can reproduce the 'H„g\ l _ gg function in [8[. Note that 



is universal, while C qq and Haa^aa contain both 



qq*-qq 



process- independent and process-dependent parts [7J. It 
is straightforward to compute the C qq and Hqq^qq func- 
tions from our results up to NNLO for any qq initiated 
process given the knowledge of the two-loop virtual cor- 
rections. 

In conclusion, we have calculated the perturbative 
quark-to-quark transverse PDF at NNLO based on a 
gauge invariant operator definition with an analytic reg- 
ulator. We demonstrate for the first time that such a 
definition works beyond the first non-trivial order. We 



extract from our calculation the coefficient functions rele- 
vant for a N 3 LL Qt resummation. Our results can be ap- 
plied to all quark-antiquark annihilation processes yield- 
ing a colorless final state, provided the NNLO virtual 
corrections are known. Combined with the recent work 
21| , our results could also be applied to the qq — > ti pro- 
cess. Our method of calculation can be easily extended 
to all parton combinations, which will be presented in a 
forthcoming article. 
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